Are there compact heavy four-quark bound states? 
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Abstract 

We present an exact method to study four-quark systems based on the hyperspherical harmonics 
formalism. We apply it to several physical systems of interest containing two heavy and two light 
quarks using different quark-quark potentials. Our conclusions mark the boundaries for the possible 
existence of compact, non-molecular, four-quark bound states. While QQnn states may be stable 
in nature, the stability of QQnn states would imply the existence of quark correlations not taken 
into account by simple quark dynamical models. 



PACS numbers: 21.45. +v, 31.15.Ja, 14.40.Lb, 12.39. Jh 



The discoveries on several fronts [1], of unusual charmonium states like X(3872) and 
F(4260) and open-charm mesons with unexpected masses like _D*j(2317) and Dq(2308), have 
re-invigorated the study of the hadron spectra. Their anomalous nature has triggered several 
interpretations, among them, the existence of four-quark states or meson-meson molecules. 
This challenging situation resembles the long-standing problem of the light-scalar mesons, 
where it has been suggested that some resonances may not be ordinary qq states, though 
there is little agreement on what they actually are 0]. In this case, four-quark states have 
been justified to coexist with qq states because they can couple to J PC = ++ without 
orbital excitation 0]. 

Any debate on the possible multiquark structure of meson resonances should be based 
on our capability to find an exact solution of the four-body problem Q]. Theoretical pre- 
dictions often differ because of the approximation method used. A powerful tool to solve a 
few-particle system is to expand the trial wave function in terms of hyperspherical harmon- 
ics (HH) basis functions. This method has been proven to be rather powerful to solve the 
nuclear [5( four-body problem. In this work we use a generalization of the HH formalism to 
study four-quark systems in an exact way. There are two basic difficulties for constructing 
HH functions of proper symmetry for a system of identical particles: first, the simultaneous 
treatment of particles and antiparticles, and second the additional color and flavor degrees 
of freedom. The method will be tested by comparing with the existing results based on 
different approximate solutions, thus establishing the validity of such approximations. Due 
to their actual interest and having in mind that systems with unequal masses are more 
promising to be bound pj, we will center on the study of QQnn and QQnn states (n stands 
for a light quark and Q for a heavy one). We will analyze the possible existence of com- 
pact four-quark bound states using two standard quark-quark interactions, a Bhaduri-like 
potential (BCN) ^ and a constituent quark model considering boson exchanges (CQC) Q. 
Both interactions give a reasonable description of the meson and the baryon spectroscopy, 
a thoughtful requirement considering that in the tetraquarks qq and qq interactions will 
contribute. 

Within the HH expansion, the four-quark wave function can be written as a sum of outer 
products of color, isospin, spin and configuration terms 

\4>cisr) = |Color)|Isospin) [|Spin> ® \R)] JM , (1) 

such that the four-quark state is a color singlet with well defined parity, isospin and total 
angular momentum. In the following we shall assume that particles 1 and 2 are the Q-quarks 
and particles 3 and 4 are the n-quarks. In the QQnn case particles 1 and 2 are identical, 
and so are 3 and 4. Consequently, the Pauli principle leads to the following conditions, 

Pi2\4>ciSR) = P-zMcisr) = -\4>cisr) , (2) 

Pij being the permutation operator of particles i and j. 

Coupling the color states of two quarks (antiquarks) can yield two possible represen- 
tations, the symmetric 6-dimensional, 6 (6), and the antisymmetric 3-dimensional, 3 (3). 
Coupling the color states of the quark pair with that of the antiquark pair must yield a color 
singlet. Thus, there are only two possible color states for a QQqq system H, 

(Color) = {|3 12 3 34 >, |6i 2 6 34 )}. (3) 

These states have well defined symmetry under permutations, Eq. (J2J). The spin states with 
such symmetry can be obtained in the following way, 

(Spin) = \((s 1 ,s 2 )S 12 ,(s 3 ,s,)S u )S) = \(S 12 Su)S) . (4) 
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The same holds for the isospin, |Isospin) = 1(^3,^4)134), which applies only to the n-quarks, 
thus I = 134. 

As mentioned above, we use the HH expansion to describe the spatial part of the wave 
function. We choose for convenience the If -type Jacobi coordinates, 

% = ^i,2(r* 2 -n) , 

fm 3 r 3 + m A r 4 m 1 ri+m 2 r 2 \ 
V2 = ^12,34 , 

V m u m 12 J 

= ^3,4(^4 - r 3 ) , (5) 

were = rrii+mj, = yfmimjjrriij, and m\ 2 ^ = mi +m2 +7723 +7724. Using these vectors, 
it is easy to obtain basis functions that have well defined symmetry under permutations of 
the pairs (12) and (34). In the HH formalism the three Jacobi vectors are transformed into 
a single length variable, p = \/ n\ + 77I + 773, and 8-angular variables, Q, that represent the 
location on the 8-dimensional sphere. The spatial basis states are given by 

( P n\R) = u n ( P )y [K] (n) , (6) 

were y\K\ are the HH functions, and [K] = {KKi 2 LMlLi 2 £ 3 £ 2 £i}. The quantum number K 
is the grand angular momentum, LMl are the usual orbital angular momentum quantum 
numbers, and li is the angular momentum associated with the Jacobi vector r/j. The quan- 
tum numbers lTi2, vL 12 correspond to the intermediate coupling of rj 1 and r) 2 . The Laguerre 
functions are used as the hyper-radial basis functions U n (p). 

The Pauli principle, Eq. (jSJ), leads to the following restrictions on the allowed combina- 
tions of basis states: 

(i) (-l) s i2+ £ i = +1, (-1)^34+1+^3 = _i f or trie |6 12 6 34 ) color state, 

(ii) (-l) s i2+^i = -1, (_i)S 3 4+/+^ 3 = +1 f or trie |3 12 3 34 > state. 

In the QQnn case the particle 2 is the antiparticle of the particle 1, and the particle 4 is the 
antiparticle of the particle 3. Assuming that C— parity is a good symmetry of QCD we can 
regard quarks and antiquarks as identical particles and impose the symmetry condition, Eq. 
(J2J), on the QQnn system as well. Coupling the color states of a quark and an antiquark can 
yield two possible representations: the singlet and the octet. These representations should 
be combined in the following way [8] {II12I34), 1812,834)} to yield a total color singlet state. 
However, these states have not definite symmetry under particle permutations (12) and (34). 
To construct symmetrized states for the QQ pair we consider the following combinations, 

\C\?) = ±=(\C l2 )+T l2 \C 21 ), (7) 

were C12 = {Ii2 ; 8i2}, and Ti 2 = +1 for a symmetric combination and —1 for an antisym- 
metric one. For light quarks the color and isospin states should be combined together to 
form states with well defined symmetry. For I z = 0, for instance, these states take the form, 

|(C 34 /34) r34 > = +\ [IC34) (\UU) ± \dd)) + r 3 4|C 4 3> (\UU) ± \dd))] , (8) 

were the plus sign stands for / 34 = state and the minus sign for the 7 3 4 = 1 state. As 
before, C 34 stands for either the singlet or the octet representations. The total color-isospin 
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states, |C , f 2 12 (C , 34 /34) r34 ) are not only good symmetry states, but also good C— parity states 
with, C = Ti2r34. Imposing the Pauli principle for the QQnn system we get the following 
restrictions: ri 2 (-l) 5l2+£l = +1, r 34 (-l) 534+ ^ = +1, on the basis states. 

Assuming non-relativistic quantum mechanics we solve the four-body Schrodinger equa- 
tion using the basis states described above. The grand angular momentum K is the main 
quantum number in our expansion and the calculation is truncated at some K value. As 
mentioned above, for our study we will use two standard quark potential models providing a 
reasonable description of the hadron spectra. In the following we draw the basic properties 
of the interacting potentials. 

The BCN model was proposed in the early 80 's by Bhaduri et al. in an attempt to 
obtain a unified description of meson and baryon spectroscopy Ifl . It was later on applied 
to study the baryon spectra 0] and four-quark (qqqq) systems [Ic|. The model retains the 
most important terms of the one-gluon exchange interaction proposed by de Rujula et al. 
TIT ] , namely coulomb and spin-spin terms, and a linear confining potential, having the form 



^ fr - K K p- r ij/ r \ 

V{r l3 ) = ■ \$ x % - ± - D + -^ C -^ r {a l ■ a,) , (9) 

16 J \ GT Tij miTTlj T ijfQ J 

where <jj are the Pauli matrices and A? are the SU(3) color matrices. The parameters 
k = 102.67 MeV fm, £>=913.5 MeV, a =0.0326 MeV" 1 / 2 fm 1 / 2 , r = 2.2 fm, m u4 =337 
MeV, and m c = 1870 MeV are taken from Ref. (To). 

The CQC model was proposed in the early 90's in an attempt to obtain a simultaneous 
description of the nucleon-nucleon interaction and the baryon spectra 12{|. It was later on 
generalized to all flavor sectors giving a reasonable description of the meson 0] and baryon 
spectra |l3| . The possible existence of four-quark states within this model has also been 
addressed [14j, Il5 



The model is based on the assumption that the light-quark constituent mass appears 
because of the spontaneous breaking of the original SU (3) i ® SU{?>)r chiral symmetry at 
some momentum scale. In this domain of momenta, quarks interact through Goldstone 
boson exchange potentials. QCD perturbative effects are taken into account through the 
one-gluon-exchange (OGE) potential as the one used in the BCN model. Finally, it incor- 
porates confinement as dictated by unquenched lattice calculations predicting, for heavy 
quarks, a screening effect on the linearly dependent interquark potential when increasing 
the interquark distance ljl . 



The model parameters have been taken from Ref. [?J with the exception of the OGE 
regularization parameter. This parameter, taking into account the size of the system, was 



fitted for four-quark states in the description of the light scalar sector [14] . being f = 0.18 
fm for mesons and f = 0.38 fm for four-quark systems. Let us also notice that the CQC 
model contains an interaction generating flavor mixing between nn and ss components. It 
allows to exactly reproduce the masses of the r\ and 7/ mesons. In the four-quark case this 
contribution would mix isospin zero QQnn and QQss components. Such contributions were 



explicitly evaluated in the variational approach of Ref. [14| for the light isocalar tetraquarks, 
giving a negligible effect. In order to make a proper comparison between thresholds and four- 
quark states we have recalculated the meson spectra of Ref. 0] with the same fo value and 
interaction used in the four-quark calculation, neglecting therefore the flavor-mixing terms. 
Explicit expressions ofthe interacting potentials and a more detailed discussion of the model 
can be found in Ref. 
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TABLE I: Energy (MeV) of L = ccnn states. 



(s,i) 


Ref. [17] HH(£j = 


0) HH 


(0,1) 


4155 


4154 


3911 


(1,0) 


3927 


3926 


3860 


(1,1) 


4176 


4175 


3975 


(2,1) 


4195 


4193 


4031 



Let us first analyze the numerical capability of the designed method to capture the 
properties of the four-quark systems. In Table [I] we present the results for different L = 
spin-isospin ccnn states calculated with the CQC model. We quote in the first column the 
results obtained within a variational calculation using Gaussian trial wave functions only 
with quadratic terms in the Jacobi coordinates [l7| • This approximation would correspond in 
our formalism to set £i = for the three Jacobi vectors. These results are given in the second 
column, reproducing exactly the variational results. The validity of this approximation can 
be judged by looking at the last column where we give the exact HH results, truncated 
at K = 24. In some cases the difference between the i± = approximation and the true 
ground state can be as large as 200 MeV. We have also reproduced the calculation of the 
(S, I) = (1,0) L = ccnn state of Refs. 0, using the BCN model. For K = 24 



we have obtained an energy of 3899.7 MeV as compared to 3904.7 MeV of Ref. [18J and 
3931.0 MeV of Ref. flof . Ref. [l8[ designed a powerful method, similar to the stochastic 
variational approach [19J, to study this particular system. Although their results are not 
fully converged, the close agreement gives confidence on both calculations. The results of 
Ref. were obtained using diagonalization in harmonic oscillator (HO) basis up to to 
N = 8. The quality of this last procedure can be judged by looking at Table [TT1 where we 



compare, for different ccnn states, the results of 10] to the HH results with K < 8 and to 
the converged HH results obtained with K restricted by our computational capability, i.e. 
K = 22 or K = 26. As can be seen the results with the basis truncated at K = 8 are similar 
to the HO results, but rather far from the converged ones. 

In spite of the shortcomings of the methods used to study four-quark systems, in the 
past, many four-quark bound states have been suggested. To analyze their stability against 
dissociation, the parity and the total angular momentum must be preserved. Additionally, 
C— parity is a good quantum number for ccnn and the Pauli principle must be fulfilled in the 
ccnn case. The thresholds can be evaluated by adding the meson masses of the dissociation 
channel. A four-quark state will be stable under strong interaction, and therefore narrow, if 
its mass lies below all allowed two-meson thresholds. Sometimes, the results of four-quark 
calculations have been directly compared to the experimental thresholds. In this case one 
could misidentify scattering wave functions as bound states. When they are referred to the 
threshold within the same model, we will see how the theoretical predictions do not imply 
an abundance of multiquark states in the data. 



TABLE II: Energy (MeV) of L = ccnn states. 



J p 


Ref. [lpj HH(K = 8) HH(K max ) 


0+ 

1+ 


3409 3380 3249 (26) 
3468 3436 3319 (22) 
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TABLE III: Energy (MeV) and probability of the different color components as a function of K 
for the ccnn J PC = 1 ++ state both for CQC and BCN models. The last rows indicate the lowest 
theoretical two-meson thresholds. The notation |g (|p) stands for relative 5-wave (P-wave). 





CQC 


BCN 


K 


E Pll Pss 


E Pn Pg8 





4141 1.0000 0.0000 


4196 1.0000 0.0000 


2 


3985 0.9822 0.0178 


4053 0.9462 0.0538 


4 


3911 0.9789 0.0211 


3994 0.9233 0.0767 


6 


3870 0.9834 0.0166 


3963 0.9236 0.0764 


8 


3845 0.9871 0.0129 


3944 0.9303 0.0697 


10 


3827 0.9905 0.0095 


3932 0.9426 0.0574 


12 


3814 0.9926 0.0074 


3920 0.9927 0.0073 


14 


3805 0.9943 0.0057 


3887 0.9990 0.0010 


16 


3797 0.9954 0.0046 


3861 0.9994 0.0006 


18 


3791 0.9962 0.0038 


3840 0.9995 0.0005 


20 


3786 0.9968 0.0032 


3822 0.9996 0.0004 


22 




3808 0.9997 0.0003 


J/ipu\s 


3745 1 


3874 1 


XcJV\p 


4281 1 


3655 1 



Let us make a brief comment on the convergence of the HH expansion. In some cases the 
convergence is slow, as can be seen by comparing Tables HLT1 and HVl This is a consequence of 
the HH formalism being better suited to treat with bound states, and most of the four-quark 
states one deals with are above the corresponding two-meson threshold, as can be seen in 
Table IIHI Due to this slowness, our calculation cannot definitively exclude large molecular 
objects (sizes of the order of 1.5-2 fm) with binding energies smaller than 1-2 MeV induced 
by long-range interactions like, for instance, the one-pion exchange tail [20]. 

Once the method has been designed, tested, and its accuracy established, we concentrate 

on a hot subject: the ccnn system as a potential structure for the X(3872). To make 

the physics clear we will compare with the ccnn system. In particular, we focus on the 
jpc = 1++ 

ccnn and J = 1 + ccnn quantum numbers to illustrate their similitude and 
differences. A complete study of all the quantum numbers will be reported elsewhere. 
The results are shown in Tables IIHI and IIVI up to the maximum value of K within our 
computational capabilities. Since we are using a complete set of HH basis, all possible 
configurations are considered in both cases. For the ccnn system, independently of the 
quark-quark interaction, the system evolves to a well separated two-meson state. This is 
clearly seen in the energy, approaching the corresponding two free-meson threshold, but also 
in the probabilities of the different color components of the wave function and in the radius. 
We denote by P n (Pss) the probability of a singlet-singlet (octet-octet) color component 
in the (cc)(nn) (or (cn)(cn)) coupling. We observe how the system evolves to two singlet 
color mesons, whose separation increases with K, dashed line in Fig. [TJ Comparing the 
theoretical predictions with the experimental threshold, Mj/^u\ 3 = 3879.57±0.13 MeV, one 
could be tempted to claim for the existence of a bound state. However, the experimental 
threshold is not reproduced by the effective Hamiltonians. In fact, in the BCN model the 
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TABLE IV: Same as Table ITU for the ccnn J p = 1+ state. 





CQC 


BCN 


K 


E 




Pu 




P&8 


E 




Pu 




-^88 





4109 





3351 





6649 


4100 





3446 





6554 


2 


3990 





3483 





6517 


3999 





3744 





6256 


4 


3931 





3577 





6423 


3954 





3981 





6019 


6 


3903 





3641 





6359 


3933 





4170 





5830 


8 


3887 





3681 





6319 


3921 





4302 





5698 


10 


3878 





3705 





6295 


3914 





4403 





5597 


12 


3872 





3720 





6280 


3910 





4478 





5522 


14 


3868 





3730 





6270 


3907 





4536 





5464 


16 


3866 





3737 





6263 


3904 





4581 





5419 


18 


3864 





3741 





6259 


3903 





4618 





5382 


20 


3862 










3901 





4647 





5353 


22 


3861 










3900 










24 


3861 










3900 










DD*\ S 


3937 




1 







3906 




1 








sum of the masses of the two mesons J/ipuo is even larger than that of XcjVi leading to a 
completely different threshold for the 1 ++ system. Thus, in any manner one can claim for 
the existence of a bound state. Similar conclusions are drawn for all quantum numbers of 
this system. 

A completely different behavior is observed in Table IIVI Here, the energy is quickly 
stabilized below the theoretical threshold. Besides, the radius is also stable, solid line in 
Fig. dj and it is smaller than the sum of the radius of the two-meson threshold. We obtain 
r 4g = 0.37 fm compared to r Ml + r M2 = 0.44 fm. We also notice a different solution for the 
probability of the color components. However, one should not directly conclude the presence 
of octet-octet components, because the octet-octet color component in the (cin 3 )(c 2 n4) basis 
can be re-expressed as a singlet-singlet color component in the (cin 4 ) (c 2 n 3 ) coupling, being 
the same physical system due to the identity of the two quarks and the two antiquarks. 
Although in the BCN model the system is slightly bound, the structure of the bound state 
is manifest for low values of K, leading one to conclude that the state could hardly be 
destroyed by small non-considered effects as could be, for example, relativity. The actual 
interest and the capability of some experiments [21j to detect double charm states makes 
this prediction a primary objective to help in the understanding of QCD dynamics. 

It is thus important to realize that a bound state should be pursued not only by looking 
at the energy, but also with a careful analysis of the radius and color probabilities. Besides, 
one should compare results within the same calculating framework, unless other effects, as 
discussed below, are considered. This detailed analysis allows us to distinguish between 
compact states and meson-meson molecules [3[ and it does consider the contribution of 



all meson- meson channels to a particular set J PC of quantum numbers 20j. Inherent to 
our discussion is a much richer decay spectrum of compact states due to the presence of 
octet-octet color components in their wave function. 

Let us notice that there is an important difference between the two physical systems 
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FIG. 1: Evolution with K of the radius (RMS) If the ccnn J PC = 1 ++ state (dashed line) and the 
ccnn J p = 1 + state (solid line) for the CQC model. 



studied. While for the ccnn there are two allowed physical decay channels, (cc)(nn) and 
(cn)(cn), for the ccnn only one physical system contains the possible final states, (cn)(cn). 
This has important consequences if both systems (two- and four-quark states) are described 
within the same two-body Hamiltonian, the ccnn will hardly present bound states, because 
the system will reorder itself to become the lightest two- meson state, either (cc)(nn) or 
(cn)(cn). In other words, if the attraction is provided by the interaction between particles 
i and j, it does also contribute to the asymptotic two-meson state. This does not happen 
for the ccnn if the interaction between, for example, the two quarks is strongly attractive. 
In this case there is no asymptotic two-meson state including such attraction, and therefore 
the system will bind. 

Therefore, our conclusions can be made more general. If we have an iV-quark system 
described by two-body interactions in such a way that there exists a subset of quarks that 
cannot make up a physical subsystem, then one may expect the existence of iV-quark bound 
states by means of central two-body potentials. If this is not true one will hardly find 
iV— quark bound states (22|. For the particular case of the tetraquarks, this conclusion is 
exact if the confinement is described by the first SU(3) Casimir operator, because when 
the system is split into two-mesons the confining contribution from the two isolated mesons 
is the same as in the four-quark system. The contribution of three-body color forces 23 
would interfere in the simple comparison of the asymptotic and the compact states. Another 
possibility in the same line would be a modification of the Hilbert space. If for some reason 
particular components of the four-quark system (diquarks) would be favored against others, 
the system could be compact j24j]. Lattice QCD calculations 25] confirm thephenomeno- 
logical expectation that QCD dynamics favors the formation of good diquarks |3|], i.e., in the 
scalar positive parity channel. However, they are large objects whose relevance to hadron 
structure is still under study. All these alternatives will allow to manage the four-quark 
system without affecting the threshold and thus they may allow to generate any solution. 



Let us finally note that in Ref. [26| the stability of the QQnn and QQnn systems has been 



analyzed in a simple string model considering only a multiquark confining interaction given 
by the minimum of a flip-flop or a butterfly potential. The ground state of systems made 
of two quarks and two antiquarks of equal masses was found to be below the dissociation 
threshold. While for the flavor exotic QQnn the binding increases when increasing the mass 
ratio mQ/m q , for the cryptoexotic QQnn the effect of symmetry breaking is the opposite, 
the system being unbound whenever mg/m, > 1. Although more realistic calculations are 
needed before establishing a definitive conclusion, the findings of Ref. 26| strengthened our 
results. 
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